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ABSTRACT 

Classical novae are important producers of radioactive nuclei, such as ^Be, ^'^N, ^*F, 
^^Na and ^^Al. The disintegration of these nuclei produces positrons (except for ^Be) 
that through annihilation with electrons produce photons of energies 511 keV and 
below. Furthermore, ^Be and ^^Na decay producing photons with energies of 478 keV 
and 1275 keV, respectively, well in the 7-ray domain. Therefore, novae are potential 
sources of 7-ray emission. We have developed two codes in order to analyze carefully 
the 7-ray emission of individual classical novae: a hydrodynamical one, which follows 
both the accretion and the explosion stages, and a Monte-Carlo one, able to treat 
both the production and the transfer of 7-ray photons. Both codes have been coupled 
in order to simulate realistic explosions. The properties of 7-ray spectra and 7-ray 
light curves (for the continuum and for the lines at 511, 478 and 1275 keV) have 
been analyzed, with a special emphasis on the difference between carbon-oxygen and 
oxygen-neon novae. Predictions of detectability of individual novae by the future SPI 
spectrometer on board the INTEGRAL satellite are made. Concerning ^^Al, its decay 
produces photons of 1809 keV but it occurs on a timescale much longer than the 
typical time interval between nova outbursts in the Galaxy, making it undetectable 
in individual novae. The accumulated emission of ^^Al from many Galactic novae has 
not been modeled in this paper. 

Key words: gamma-rays: general - novae, cataclysmic variables - nucleosynthesis - 
stars: abundances - white dwarfs 
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1 INTRODUCTION 

Thermonuclear runaways on accreting white dwarfs are at 
the origin of nova explosions. Some unstable proton-rich 
nuclei, such as ^^N, ^*0, ^^O, ^'^F and ^*F, are produced 
during hydrogen-burning. These nuclei are /3+-unstable and 
therefore decay through the emission of positrons. This has 
two important consequences. First of all, the decay times 
of "O (r=102s), ^^O (r=176s) and ^''F (r=93s) are such 
that these nuclei can be transported by convection to the 
outer layers of the envelope (during the runaway) before 
decaying; their subsequent decay in these outer shells of 
the envelope is largely responsible for the expansion and 
of the large increase in luminosity of the nova. Second, the 
positrons emitted by all those nuclei annihilate with elec- 
trons leading to the emission of photons with energy equal 
or below 511 keV. (r=862s) and ^**F (r=158min) are 
the most important contributors to 7-ray emission, since 
their decay timescales are larger and they decay when the 
envelope is transparent enough. Other radioactive elements 
are synthesized during nova explosions: carbon-oxygen (CO) 
novae are important producers of ^Be (T=77days) (see Her- 
nanz et al. 1996), which emits a photon of 478 keV when it 



decays through an electron capture, whereas oxygen-neon- 
magnesium novae (ONeMg) produce ^^Na (r=3.75yr) and 
^''Al (r=1.04x lO'^yr), which decay emitting a positron and 
a photon of 1275 and 1809 keV, respectively. 

Some previous works have pointed out the potential im- 
portance of classical novae as 7-ray emitters (Clayton 1981, 
Clayton & Hoyle 1974, Leising & Clayton 1987) . Concerning 
the nucleosynthesis of radioactive nuclei in novae, it has been 
analyzed by means of parametrized studies (Weiss & Tru- 
ran 1990, Nofar, Shaviv & Starrfield 1991), and by complete 
hydrodynamic simulations (see Politano et al. 1995 and Pri- 
alnik & Kovetz 1997, for recent works concerning ONeMg 
and CO novae, respectively, and also Jose & Hernanz 1997). 
But there is a lack of detailed studies of the 7-ray output of 
novae and its relation with the production of radioactive nu- 
clei. Only the use of realistic profiles of velocities, densities 
and chemical abundances, necessary to know the production 
and transfer of 7-rays during nova explosions, will provide 
the correct 7-ray spectra of classical novae (Hernanz et al. 
1997a and b). 

From the observational point of view, two types of ob- 
servations have to be distinguished: individual and cumu- 
lative. The emission of short and medium-lived radioactive 
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nuclei (478, 511 and 1275 keV lines and continuum) can be 
detected in individual novae, provided they are close enough. 
The emission of very long-lived ^^Al nuclei can contribute to 
the diffuse emission in the Galaxy, but can not be detected 
in particular novae. The emission of the medium-lived nuclei 
('^Be and ^'^Na at 478 and 1275 keV) can also be accumu- 
lated from different novae, as their lifetimes are longer than 
the typical time interval between two consecutive nova ex- 
plosions in the Galaxy. Up to now, no positive detection of 
novae in 7-rays has been produced, neither for particular 
novae in the 478 or 1275 keV lines (Harris, Leising & Share 
1991, lyudin et al. 1995), nor for the cumulative emission of 
novae at the Galactic center in these same lines (Leising et 
al. 1988, Harris et al. 1991 and 1996). 

We will concentrate in this paper on the 7-ray emission 
from individual novae, by means of hydrodynamical simula- 
tions which include both the accretion and explosion stages, 
with complete nucleosynthesis, and a Monte-Carlo technique 
to follow the production and transfer of 7-ray photons. The 
analysis of the cumulative emission from novae in the Galaxy 
is currently under way. The nova models are presented in 
section 2, with special emphasis on the synthesis of radioac- 
tive nuclei and their density and velocity profiles. In section 
3 we describe the Monte-Carlo simulation technique and in 
section 4 its results, in the form of 7-ray spectra and light 
curves for the continuum and the main lines of interest. Fi- 
nally, our discussion and conclusions follow. 



2 NUCLEOSYNTHESIS OF RADIOACTIVE 
NUCLEI 

The accretion and explosion stages have been followed by 
means of an implicit hydrodynamic code, which includes 
a complete nuclear reaction network; one hundred nuclei, 
ranging from to ''"Ca, linked through 370 nuclear reac- 
tions, have been included (see Jose 1996 and Jose & Hernanz 
1997a and b, for details). A time-dependent formalism for 
convective transport has been included whenever the char- 
acteristic convective timescale becomes larger than the in- 
tegration time step (Wood 1974). Partial mixing between 
adjacent shells is treated as a diffusion process (Prialnik, 
Shara & Shaviv 1979). CO as weU as ONeMg novae have 
been computed, with an accretion rate of 2x10"^" MQyr~^ 
(other values for the accretion rate have also been consid- 
ered, but they do not introduce noticeable changes in the 
properties of the radioactive nuclei synthesized). The initial 
composition of nova envelopes is far from being understood, 
but some enrichment with material from the white dwarf 
core is needed in order to explain both the explosion mech- 
anism itself and some of the observed abundances in novae. 
Thus, as the real mechanism of enhancement is not clear at 
present, a compromise solution consisting of a parametrized 
mixture between the outermost shells of the underlying core 
and the solar-like accreted matter is adopted. We have as- 
sumed a 50% degree of mixing, although other possibilities 
have also been analyzed. For the CO novae, the adopted 
composition of the underlying core consists of a 50%-50% 
mixture of ^^C and of ^^O, whereas for ONeMg novae (where 
the initial composition is more crucial for the synthesis of 
the radioactive nucleus ^^Na) we have adopted a detailed 
outer core initial composition from Ritossa, Garcfa-Berro 



& Iben (1996) (indicating that these white dwarfs are al- 
most devoid of magnesium, which makes it more correct to 
call them oxygen-neon (ONe) white dwarfs). The exact ini- 
tial composition is an important issue, since the presence of 
■^"Ne and its amount, as well as that of magnesium, is cru- 
cial for the final amount of the radioactive nuclei ^'^Na and 
^^Al. In particular, the new initial abundances we are using 
favor a lower ^^Na and ^^Al production than the ones com- 
monly adopted, based on old nucleosynthesis calculations by 
Arnett & Truran (1969) (Jose, Hernanz & Coc 1997, Jose 
& Hernanz 1997a and b). Recent observations of individual 
novae with the COMPTEL instrument on board the COMP- 
TON Gamma-Ray Observatory (CGRO) have not detected 
the 1275 keV emission in any of the observed novae. Only 
upper limits to the flux, which translate into upper lim- 
its to the ^^Na ejected mass, were obtained (lyudin et al. 
1995). COMPTEL observations have also provided detailed 
maps of the Galactic 1809 keV emission, corresponding to 
^^Al decay (Diehl et al. 1995). The main features of these 
maps seem to indicate that this emission is rather correlated 
with a young population, thus discarding novae as the main 
contributors to it (see Prantzos & Diehl 1996 for a recent 
review). Our low ^^Na and ^®A1 yields are consistent with 
both observational results. 

The production of radioactive nuclei, as well as some 
other relevant properties, of some selected CO and ONe no- 
vae are displayed in Table ^ (see Jose and Hernanz 1997a 
and b for more details). Models Al and A2 correspond to 
ONe novae with masses 1.15 Mq and 1.25 Mq, respectively, 
whereas models Bl and B2 correspond to CO novae, with 
masses 0.8 and 1.15 Mq, respectively. Models Al, A2 and 
B2 share some general properties: they have similar ejected 
masses (~ 10~^ Mq) and similar kinetic energies (~ 3 x 10^^ 
erg.g~^), with typical average velocities between 2000 and 
3000 km.s'^. Concerning model Bl, its lower initial mass 
leads to a higher mass of the accreted envelope and thus to 
a higher ejected mass but with a lower mean velocity (~ 
1000 km.s~^), which translates into a lower kinetic energy. 
This makes this model more opaque to the 7-rays, as will be 
shown below. Concerning the synthesis of radioactive nuclei, 
the most important result is that ^Be is mainly produced in 
CO novae whereas ^^Na is only produced in ONe novae. This 
result is a consequence of the different initial compositions of 
CO and ONe novae, as discussed above. The complete pro- 
files of density, velocity and chemical composition obtained 
with the hydrodynamical code, taken when the unacceler- 
ated expansion begins, have been used as inputs for the cal- 
culations performed with the Monte-Carlo code developed 
to obtain the gamma-ray emission of those classical novae. 



3 GAMMA-RAY PRODUCTION AND 
PROPAGATION: MONTE-CARLO 
SIMULATION 

In order to simulate the evolution of the 7-ray spectra of 
the different models of novae considered, we have developed 
a 7-ray transfer code on the basis of the method described 
by Pozdnyakov, Sobol & Sunyaev (1983) and Ambwani & 
Sutherland (1988). This code has been successfully applied 
to the study of the spectra of type la supernovae (Gomez- 
Gomar, Isern & Jean 1997). It is based on the Monte-Carlo 
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Table 1. Main properties of the ejecta one hour after peak temperature. 



Model 


Mwd 


Mejcc 


<Ek > 


^Be 




13 N 




18p 


22 Na 


Al 


1.15 


1.8 10-5 


3.1 10^^ 


~ 




5.5 10- 


-9 


7.1 10-* 


9.8 10-10 


A2 


1.25 


1.6 10-5 


3.3 lO^^ 


1.2 10- 


ll 


2.9 10- 


'8 


6.7 10-« 


1.6 10-9 


Bl 


0.8 


6.3 10-5 


8 10^5 


7.8 10- 


■ 11 


1.6 10" 


-7 


1.7 10-^ 


~ 


B2 


1.15 


1.4 10-5 


3.2 lO^"^ 


1.1 lo- 


■10 


1.3 10" 


-8 


3.6 10-* 


~ 



M„di total ejected mass (Mcjcc) and ejected mass of each isotope are in Mq, 
and the average kinetic energy < Ek > is in erg g-^. Models Al and A2 
correspond to ONe novae, whereas Bl and B2 are CO ones. 



technique, which allows to treat the comptonization of high 
energy photons without approximations. In the code, the 
radioactive decays of ^Be, ^^N, ^^F and 22Na are included 
to generate the initial photons. A comment must be made 
regarding '^Be. This nucleus decays via the capture of an 
electron of the K-shell. This implies that at high tempera- 
tures its decay time in the laboratory could change because 
of ionization. However, as the lifetime of this nucleus is quite 
long we assume that for most of the time the temperature of 
the envelope is low enough to avoid complete ionization and 
the decay time is always taken to be equal to the laboratory 
time. 

Once photons are generated according to the relative 
isotopic abundances and rates of disintegration of the above 
mentioned radioactive nuclei, their trip across the expanding 
ejecta is simulated by taking into account the three differ- 
ent interactions which affect their propagation, i.e. Comp- 
ton scattering, photoelectric absorption and production of 
e~—e'^ pairs. The cross sections for all these interactions is 
computed taking into account the precise composition of the 
ejecta. In the case of Compton scattering the cross section 
adopted is the usual Klein-Nishina expression, while absorp- 
tion and pair production cross sections are taken from the 
compilation of experimentally evaluated data maintained by 
the Brookhaven National Laboratory. In the code the effects 
of the expansion velocity are included an thus realistic line 
profiles are obtained. Besides direct emission by radioactive 
decay, photons are also emitted by positron annihilation. 

The treatment of positron annihilation has required 
particular attention, since during some phases of the nova 
evolution this is the dominant emission mechanism, being 
"N, and ^^Na the main contributors. When a positron 
is emitted, either it escapes without interaction with the 
ejecta or it annihilates. We have considered that in a nova 
envelope positrons thermalize before annihilating. This ap- 
proximation is wrong in less than 1% of the cases in an 
electronic plasma (Leising & Clayton 1987). Concerning a 
neutral envelope, the excitation cross section dominates any 
other interaction above energies of ~100 eV (Bussard, Ra- 
maty & Drachman 1979) and thus positrons lose energy un- 
til they reach this extremely low value. To reproduce this 
braking effect in the simulations, positrons propagate until 
they cross an equivalent column mass density of 0.2 g cm-2 
(measured in straight line) (Chan & Lingenfelter 1993). This 
is the mean range expected for a 0.6 MeV positron braked 
down to energies ~100 eV through elastic collisions with the 
surrounding medium, and neglecting the effects of the mag- 
netic field on its propagation. When a positron is thermal- 
ized, it covers a negligible distance and then annihilates. For 
temperatures and densities of a nova envelope the positrons 
form positronium in ~90 % of annihilations (Leising & Clay- 



ton 1987) while in the remaining 10 % of cases they annihi- 
late directly. 25% of times positronium is formed in singlet 
state, which leads to the emission of two 511 keV photons, 
while 75% of times it is formed in triplet state, leading to 
a three photon annihilation. The continuum spectrum of 
photons produced by triplet annihilations is obtained from 
the results of Ore & Powell (1949). To summarize, once a 
positron is produced we follow its trip until it escapes or 
covers the average energy loss distance. In the latter case 
it produces positronium 90% of times leading to triplet to 
singlet annihilations in the proportion 3:1, while in 10 % of 
the cases it annihilates directly. 



4 GAMMA-RAY SPECTRA AND LIGHT 
CURVES 

The 7-ray signatures of CO and ONe novae are very differ- 
ent. In the CO case, an important production of ^Be ensues 
(Hernanz et al. 1996), leading to line emission at 478 keV. 
There is also line emission at 511 keV, related to e--e^ 
annihilation, and a continuum related to comptonized 511 
keV emission and positronium decay (see below for details) . 
On the other hand, ONe novae are '^^Na producers, leading 
to line emission at 1275 keV (plus the 511 keV and below 
emission similar to the one of the CO cases). 

In figures 1, 2, 3 and 4 the spectral evolution of models 
Al, A2, Bl and B2, representing the emission expected for 
these novae at a distance of 1 kpc, is shown. These sequences 
correspond to the early stages of the nova outburst up to 2 
days after peak temperature. In all cases a continuum com- 
ponent, basically below 511 keV, as well as some lines (478 
keV, 511 keV and 1275 keV) appear. The relative intensity 
of the continuum and the lines, the particular lines appear- 
ing and the temporal evolution of the spectra depend on the 
model considered. In what follows we will discuss in detail 
the properties of the continuum and the lines. 

4.1 Continuum 

In all models a strong continuum dominates the spectrum 
during the early epochs of the expansion, when the ejecta is 
optically thick. It is produced both by the comptonization 
of photons emitted in the 511 keV line and by the photons 
directly emitted in the continuum by positronium annihila- 
tion. In order to display the evolution of the continuum in 
all models, we have computed the 7-ray fight curves for the 
emission in a band ranging from 170 keV to 470 keV (see fig. 
^ . This band contains most of the continuum and no appre- 
ciable contribution of the 478 and 511 keV lines. Further- 
more, it avoids energies where the instrumental background 
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Figure 1. Evolution of the 7-ray spectrum of model Al (1.15 
Mq, ONe nova; d=l kpc). Spectra correspond to 6, 12, 24 and 
48 hours after peak temperature. Notice that the 1275 keV line 
has not yet reached its maximum flux at t=48 h (see Table 3 and 
Figure 6). 



Figure 3. Evolution of the 7-ray spectrum of model Bl (0.8 , 
CO nova; d=l kpc). Spectra correspond to 6, 12, 24 and 48 hours 
after peak temperature. 
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Figure 2. Evolution of the 7-ray spectrum of model A2 (1.25 
Mq, ONe nova; d=l kpc). Spectra correspond to 6, 12, 24 and 
48 hours after peak temperature. 



for the SPI spectrometer on board the future INTEGRAL 
satellite will be higher (Jean 1995). The duration of this 
period is determined by the mass of the ejecta and its ex- 
pansion rate, being longer for more massive envelopes. As 
shown in figure |^, a reduction of the flux by a factor of 100 
occurs at around day 1. In this phase, the main source of 7- 
ray photons is the annihilation of positrons coming from the 
decay of ^^N and ^*F. In all cases the continuum displays a 
cutoff at low energies (see figs. 1 to 4). During the first hours, 
the cutoff is caused by photoelectric absorption, which has 
larger cross sections than Compton scattering at low ener- 
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Figure 4. Evolution of the 7-ray spectrum of model B2 (1.15 
Mq, CO nova; d=l kpc). Spectra correspond to 6, 12, 24 and 48 
hours after peak temperature. 



gies and acts as a sink of comptonized photons. Because of 
the strong effect of the composition on the cross section of 
photoelectric absorption, the cutoff is located at higher en- 
ergies for the ONe rich envelopes (~30 keV to be compared 
to ~20 keV for CO models). As the ejecta expands, the 
optical thickness of the envelope decreases and comptoniza- 
tion becomes unable to "pump" photons to low energies, 
leading to the shift of the cutoff towards higher energies 
and to smoother slopes. When the ejecta becomes transpar- 
ent, the contribution of comptonization to the continuum 
completely disappears and the continuum corresponds ex- 
clusively to the photons coming from positronium annihi- 
lation, which displays its intrinsic emission spectrum (Ore 
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Figure 5. Evolution of the 7-ray emission in the band 170 keV- 
470 kcV (d=l kpc) for the different models. Solid line corresponds 
to model Al (1.15 Mq, ONe nova), dashed line to model A2 
(1.25 M0, ONe), dotted line to model Bl (0.8 Mq, CO), and 
dot-dashed hnc to model B2 (1.15 Mq, CO). 



& Powell 1949). In the transition from an optically thick 
to an optically thin envelope, a feature appears at 170 keV 
(the energy of a backscattered 511 keV photon) which is re- 
lated to the fact that very few photons are scattered more 
than once, leading to a steep decrease of the continuum in- 
tensity at this energy (see figs. 1 to 4). In the case of CO 
novae (models Bl and B2), after the decay of ^^N and 
no relevant e^ emitters are left and the continuum prac- 
tically disappears (flux F< 10"'' photons/s/cm^) between 
days 1 and 2. However, for ONe novae (models Al and A2) 
■^^Na is a secondary source of e^, with a longer lifetime than 

and "F. It emits a e+ in 90 % of its decays. For these 
models, the continuum would in principle decline with the 
decay time of ^^Na (3.75 years); however, when the ejecta 
is so transparent that even the positrons generated by the 
decay of '^^Na can escape without interacting with matter, 
the continuum disappears. This happens at t~7.5 days for 
models Al and A2 (see flg. ^ . 

The maximum for the emission in the continuum co- 
incides with that for the 511 keV line in all cases, as we 
will see below (fig. It appears around 5-6 hours after 
peak temperature. This maximum is mainly associated with 
the decay of ^*F (r=2.6 hours). However, in models Bl and 
B2 a previous and stronger maximum appears around one 
hour after peak temperature, caused by the emission of the 
short-lived ^^N. This maximum is related to the high abun- 
dance of ^'^N in the outer layers of these models. In Table 
^ the positions and intensities of these maxima, as well as 
the moment when the continuum intensity falls below 10~® 
photons/s/cm^ (toff) are summarized. 

Although the continuum duration is very short as com- 
pared to the nominal observation times of INTEGRAL'S 
SPI instrument, its large intensity during this short pe- 
riod makes it interesting to compute the detectability dis- 
tance for the different models. This is the value appearing 
in the last column of Table and it has been computed 



only 10 hours, in order to obtain good sensitivities when the 
flux in the continuum is still high. All the observations start 
5 hours after peak temperature. All the computed models, 
except the low-mass CO nova (Bl), could be detected in the 
continuum up to the Galactic center. 



4.2 Lines 

Lines are the best candidates for the detection of novae. This 
is particularly true in the case of a high resolution spec- 
trometer like integral's SPI because a flux emitted in a 
narrow line is detected with these instruments with signifi- 
cances much higher than if it was spread in a broad band. 
The strongest 7-ray lines emitted by CO and ONe novae 
are shown in figures 1 to 4, where the spectral evolution of 
models Al, A2, Bl and B2 during the first 2 days is dis- 
played. Two phases can be distinguished in the evolution of 
all lines. During the early nova expansion (up to a week, de- 
pending on the model), the envelope is optically thick and 
the intensity is not only determined by the instantaneous 
abundance of the parent isotopes but also by their distribu- 
tion as a function of the optical depth (which depends on 
densities and expansion velocities). At this stage, the evolu- 
tion of lines is affected by two timescales: the decay times of 
the radioactive isotopes and the characteristic timescale of 
envelope expansion. This phase includes both the rise of the 
intensity of the long-lived 1275 keV and 478 keV lines, and 
the rise and decay of the 511 keV line. Later on, when the 
effects of absorption and comptonization become negligible, 
the intensity of the lines is determined by the total mass 
of the radioactive isotopes and the evolutionary timescale 
of the 1275 keV and the 478 keV lines simply corresponds 
to the decay times of ^^Na and '^Be (see figures |^ and ^ 
where light curves for the 1275 and 478 keV lines are shown 
for ONe and CO novae, respectively). However, even at late 
times line profiles (which have to be taken into account to 
determine the detectability of the lines) are determined by 
the distribution of ^^Na and ^Be. 



4.2.1 The 478 keV and 1275 keV lines 

The observation of any of these lines would allow to dis- 
criminate between ONe (Al and A2) and CO novae (Bl 
and B2). While the 1275 keV line can only be observed in 
the case of ONe novae (a negligible ^^Na amount is pro- 
duced by CO novae), the reverse is true for the 478 keV 
line, except for model A2, where a small amount of ^Be (see 
Table 1) leads to a weak line at 478 keV (see figures 1 to 4). 
In Table ^ we summarize the epoch of the maximum, the 
total line flux at maximum and the detection distances for 
the lines at 3 cr significance. These detection distances have 
been computed for INTEGRAL'S SPI assuming an obser- 
vation time of lO^s starting at maximum. In all cases, the 
detection distances are ~0.5 kpc. For the computation of the 
detection distances, we have taken into account the profiles 
of the lines: typical widths (FWHM) when the envelope is 
transparent are in the range between 3 keV (model Bl) and 
7 keV (model B2), for the 478 keV line, and around 20 keV 
for the 1275 keV one (for both Al and A2 models). 
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Table 2. Positions and intensities of the flux maxima in the continuum (170 keV-470 keV) 
for the different models. Subscript 1 corresponds to values (if they exist) for the maxima 
associated with ^^TSS decay, while subscript 2 corresponds to maxima associated with 
decay. At toff, the flux in the continuum falls below 10~® photons/s/cm-^ . D is the maximum 
detectability distance for INTEGRAL'S SPI (see text). 



Model 


ti [hours] 


Fluxi [photons/s/cm^] 


t2 [hours] 


Flux2 [photons/s/cm'^] 


toff [days] 


D [kpc] 


Al 






5.5 


7.4 10-2 


7.5 


8.8 


A2 






5 


9.1 10~2 


7.5 


8.8 


Bl 


1 


1.7 10-3 


8.4 


8.3 10-* 


1.4 


1.1 


B2 


1 


1.7 10-1 


6.2 


5.7 10-2 


1.7 


8 



The rise phase for the 1275 keV line lasts ~7 days. Soon 
after the maximum, the light curve reaches the stable decline 
phase with the expected decay time of 3.75 years. During 
this phase, the line intensities are directly determined from 
the ejected '^^Na mass of each model, which grows with the 
mass of the progenitor white dwarf (Table hi), a tendency 
that had been previously pointed out by other authors (Star- 
rfield et al. 1993, Politano et al. 1995). 

On the basis of the calculations by Starrfield et al. 
(1978), Clayton (1981) first suggested that the 478 keV line 
produced by the decay of '^Be could be detectable for a 
nearby nova. This possibility was subsequently discarded by 
the results obtained by BofBn et al. (1993), who considered 
an updated reaction network, but did one-zone model nucle- 
osynthesis calculations. The new calculations performed by 
Hernanz et al. (1996) show that convection plays an impor- 
tant role and that an appreciable amount of '^Be is produced 
by a CO nova. The possible detection of the corresponding 
478 keV line would be an important confirmation of the 
thermonuclear runaway model for novae and its associated 
nucleosynthesis. In fig. |^ we distinguish two phases in the 
evolution of the 478 keV line for the two CO models. During 
the first ^^1.5 days after peak temperature, the intensity of 
the line is completely dominated by the contribution of the 
very strong continuum generated by ^^^N and ^^^F emission. 
During this epoch, the line reaches its absolute maximum. 
Later on, the contribution of the continuum disappears and 
the line follows a "typical" light curve. The properties of 
the line at its late maximum are shown in Table ^. This late 
maximum is reached at day 5 for the massive CO nova B2 
(F=2.6 10-*^ photons/s/cm^). The low-mass CO nova, Bl, 
has a less intense maximum (F=1.4 10^'' photons/s/cm^), 
related to its lower ^Be content, appearing later (t'~13 days), 
because it is more opaque. 



4.2.2 The 511 keV line 

This line is expected to be, by far, the strongest 7-ray fea- 
ture emitted by any nova and it appears with high intensity 
in all models. However, it is much fainter in the low-mass CO 
model (Bl), which has a more opaque envelope because of its 
higher mass and lower velocity. Several isotopes contribute 
to the 511 keV emission (see section 1), with decay times 
ranging from minutes to years. This results in a relatively 
complex temporal evolution. Besides, as the line reaches its 
maximum during the 10 first hours after peak temperature, 
when the envelope is optically thick, its intensity is strongly 
affected by the details of the distribution of the radioactive 
nuclei and the expansion velocities. Thus its detection would 
provide critical information about the properties of novae 
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Figure 6. Light curves for the 1275 keV line (d=l kpc). Solid 
line corresponds to model Al (1.15 Mq, ONe nova) and dashed 
line corresponds to A2 (1.25 Mq, ONe). 
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Figure 7. Light curves for the 478 keV line (d=l kpc). Dotted 
line corresponds to model Bl (0.8 Mq, CO nova), whereas dot- 
dashed line corresponds to model B2 (1.15 Mq, CO). 
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Table 3. Properties of the 478 keV and 1275 keV lines: time and flux at maximum for an assumed distance of 1 kpc. D is the maximum 
detectabiUty distance for INTEGRAL'S SPI (see text). 



Model 


line 


t [days] 


Flux [photons/s/cm"^] 


D [kpc] 


Al 


1275 keV 


7.5 


3.8 lO-*^ 


0.5 


A2 


1275 keV 


6.5 


6.1 10"^ 


0.6 


Bl 


478 keV 


13 


1.4 lO"*' 


0.4 


B2 


478 keV 


5 


2.6 10-*^ 


0.5 



Table 4. Properties of 511 keV line: time and flux at maximum for an assumed distance of 1 kpc. D is the maximum detectability 
distance for INTEGRAL'S SPI (see text). 

Model t [hours] Flux [photons/s/cm^] D [kpc] 

Al 6 1.3 10-2 13.5 

A2 5 1.6 10-2 12.5 

Bl 1/7.5 1.7 10-3/1.4 lO-'' 1.6 

B2 1/6.5 10-2/7.3 10-3 9 
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Figure 8. Light curve evolution for the 511 keV line (d=l kpc). 
Solid line corresponds to model Al (1.15 Mq, ONe nova), dashed 
line to model A2 (1.25 Mq, ONe), dotted line to model Bl (0.8 
Mq, CO) and dot-dashed line to model B2 (1.15 Mq, CO). 



ejecta. The properties of the line in the different models ap- 
pear in Table W, whereas the light curves are displayed in 
figure ^ As happens with the continuum, the 511 keV line 
shows two maxima in models Bl and B2, which are related 
to the decay of and ^®F, while the first maximum does 
not appear in models Al and A2. The maximum intensities 
are shown in Table ^. Light curves for models Al and A2 are 
very similar, specially in the period between 0.5 days and 1 
day after peak temperature. The intensity of the late max- 
imum for model B2 is approximately 2 times smaller than 
that of models Al and A2, although later than 12 hours af- 
ter peak temperature its hght curve converges with those of 
models Al and A2. Model Bl emits lower fluxes at the two 
maxima (around 100 times smaller than ONe models for the 
second maximum) , because of the larger optical thickness at 
these early epochs. Later than 2 days after peak tempera- 
ture, models Al and A2 present a more stable emission at 
511 keV produced by the contribution of ^-^Na which emits 
a e^ in 90 % of its decays. As happens for the continuum 



(see above), this remaining emission totally disappears later 
than a week after peak temperature, when positrons emitted 
by 22Na freely escape without annihilating. 

For the estimation of the maximum distances for de- 
tectability (see Table we have adopted in this case the 
same observation parameters as for the continuum, because 
of the shortness of the 511 keV emission (for comparison, 
the sensitivity offered by OSSE, on board CGRO satellite, 
at 511 keV is approximately 10 times worse than that offered 
by INTEGRAL'S SPI for the same conditions). For models 
Al, A2 and B2, the line would be observed at distances 
equivalent to the Galactic center (13.5 kpc, 12.5 kpc and 9 
kpc) although the Galactic center itself should be avoided 
because of its 511 keV characteristic emission. The lower in- 
tensity of the 511 keV line for the low-mass CO nova, Bl, 
makes it only detectable for a distance smaller than 1.6 kpc. 
The typical width (FWHM) of the 511 keV hne is 8 keV 
(except for model Bl, which has a width of 3 keV). 

The fundamental difficulty posed by the detection of 
the 511 keV line is, of course, its short duration. The de- 
tectability of the line critically depends on the time at which 
the observation starts. To visualize this dependence, we have 
computed the maximum distance of detection for a 10 hours 
observation as a function of the starting time of the obser- 
vation (Fig. It is easy to appreciate from the figure how 
fast the chances of detection decline. The maximum distance 
of detection falls below 1 kpc if the observation starts later 
than ~20 hours after peak temperature. Thus a strategy for 
early detection of novae is necessary in order to observe this 
emission. A chance is offered by the fact that INTEGRAL 
will spend an important part of its observation time per- 
forming a survey of the Galactic plane, where most of the 
novae are expected to occur. Because of its large field of 
view (~ 16°), the possibility of an early nova detection by 
integral's SPI during the survey is not negligible. 



5 DISCUSSION AND CONCLUSIONS 

It is important to mention that all theoretical models, in- 
cluding ours, do not fit some of the ejected masses observed 
in novae eruptions (although not all observations of the same 
nova give the same value of the ejected mass). As an ex- 
ample, our ONe model Al (Mwd=l-15) fits quite well the 
observed abundances of the neon nova QU Vul 1984, but 
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Figure 9. Maximum distance of detection for the 511 keV line 
emission as a function of the starting observation time. Solid line 
corresponds to model Al (1.15 Mq, ONe nova), dashed line to 
model A2 (1.25 Mq, ONe), dotted line to model Bl (0.8 Mq, 
CO) and dot-dashed line to model B2 (1.15 Mq, CO). The results 
correspond to a 10 hours observation with INTEGRAL'S SPI. 



the ejected mass of this nova deduced from observations is 
between 2x10"^ and ~ 10"^ Mo(Saizar et al. 1992), almost 
two orders of magnitude larger than the theoretical one in 
the most unfavourable case. Thus, the flux of the 1275 keV 
line could be considerable larger for some particular ONe 
nova, but no present theoretical models are able to produce 
simultaneously such large ejected masses and neon in the 
ejecta (see Shara 1994, Shara & Prialnik 1994 and Starrfield 
et al. 1997 for recent discussions and some suggestions rela- 
tive to this problem). In view of this still open problem, we 
have made some estimations of the effect of the mass of the 
envelope, as well as the mean velocity of the ejecta, in the 
flux of the predicted 7-ray lines. Concerning the 511 keV 
line, an increase by a factor of ~20 in the ejected mass only 
reduces the flux at maximum by a factor of ~3. On the other 
hand, a reduction of the speed by a factor of ~8 reduces the 
flux by a factor of ~100. Concerning the medium and long 
lived-emission from ^Be and ^^Na (at 478 and 1275 keV, re- 
spectively), ejected mass influences in the expected way the 
fluxes (a factor of 10 increase in mass leads to a factor 10 
increase in flux), because that emission is mainly produced 
when the envelope is transparent. By the same reason, the 
velocity only influences the rise to maximum, because this 
is the phase where transparency plays an important role. 

One of the main conclusions of our study is that the 
7-ray signatures of CO and ONe novae are very different, 
except for the early emission. During the first hours after 
explosion (i.e. after peak temperature or roughly after peak 
bolometric luminosity), continuum emission below 511 keV 
dominates, as well as an intense line at 511 keV, for all the 
models. But this emission can more easily be detected "a 
posteriori", with a retrospective analysis of the data from 
instruments like BATSE on board CGRO (Fishman et al. 
1991). The reason is that this early emission happens most 
probably before the maximum in visual magnitude and even 



before detection for the majority of novae: around 2 days be- 
fore visual maximum for a very fast nova, between 2 and 5 
days for a fast nova and between 5 and 10 days before maxi- 
mum for a moderately fast nova; the majority of ONe novae 
have been classified as fast novae. We are mainly concerned 
by moderately fast to very fast novae, because these are the 
most luminous ones and thus also the most easily detectable 
either in visual or in 7-rays. Concerning the late emission, 
it is very different in a CO nova than in an ONe one. In 
the first case, emission from ^Be at 478 keV is obtained, 
with fiuxes lasting for some months but just in the limit of 
detectability by present and future instruments (see how- 
ever our comment concerning ejected masses, which applies 
both to CO and to ONe novae). A similar situation concern- 
ing detectability occurs in the ONe case, but related to the 
1275 keV line, which lasts around 3 years. We understand 
why neon novae have not been detected by COMPTEL, and 
we are moderately optimistic concerning the future SPI in- 
strument on board INTEGRAL, specially if a nearby and 
high envelope mass nova explodes during the duration of 
the mission. 
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